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Abstract 

From an analysis of the ionisation energy loss of charged particles selected 
from a sample of 147926 e"'"e~— >T"'"r~ candidates recorded in the OPAL detec- 
tor at e^e~ centre-of-mass energies near the resonance, we determine the 
branching ratios: 

Br(T" ^ K"7r"7r+(7r°)) = 0.343 ± 0.073 ± 0.031 % 
Br(T^ ^ Vr K"7r"K+(7r°)) = 0.159 ± 0.053 ± 0.020 %, 

where the (vr^) notation refers to decay modes with or without an accompanying 
7r°. The t~— ^-z^,- K~7r 7r^('?i"'^) final states occurring through r~— > i^,- K^Kg(7r^) 
are treated as background in this analysis. 

We also examine the resonant structure of t~ ^ Vt^^ candidates. Un- 

der the assumption that the resonant structure is dominated by the Ki reso- 
nances, we determine: 

Br(T- z^^Ki(1270)) 
^ - Br(r-^ z..Ki(1400))+Br(r-- z..Ki(1270)) 

In all results, the first uncertainties are statistical and the second are systematic. 
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1 Introduction 



The study of rare decays of the tau lepton has been made possible by the availabihty 
of large, low background tau-pair samples, such as those produced at the resonance 
at LEP I. This paper reports on an analysis of three-prong tau decay modes contain- 



3 



ing charged kaons, using the complete set of data collected by the OPAL experiment 
between 1990 and 1995 at e~^e~ centre-of-mass energies near the T? resonance. The 
excellent charged particle identification capability of the OPAL detector is exploited 
to study these states and obtain measurements of their branching ratios. 

The lepton can only decay to hadronic final states of strangeness or ±1. These 
two allowed values of the final state strangeness imply that if a first-order weak decay 
of a tau to a three-prong final state includes exactly one kaon, then that kaon, if it is 
charged, must have the same charge as the initial tau. Likewise, if a first-order weak 
decay of a tau to a three-prong final state includes exactly two charged kaons, those 
kaons must be of opposite charge. Thus t~^Vt K~7r^7r+ and t~^Vt K~7r~K"'" decays 
can both occur as first-order weak interactions. Decays which violate the strangeness 
constraint, such as r~— z/,- 7r~7r~K"'" and r~— ^-i/,- K~K~7r"'", can only occur as second 
order weak interactions, and are therefore highly suppressed. The charge conjugate 
decays are implied in these interactions and throughout this paper. 

Table lists three-prong decays of the tau with charged kaons in the final state, 
along with the current experimental measurements and theoretical expectations for 
each decay branching ratio [|ll-0. Due to the relative rarity of tau three-prong fi- 
nal states with charged kaons, and the experimental difficulties involved in charged 
particle identification, the branching ratios and resonant structure of these decays 
are not well known. We determine the branching ratios for r~— K~7r~7r"'"(7r'') and 
r~— > z/t- K~7r~K"'"(7r°), and, to verify that the charged particle identification techniques 
used in this analysis yield unbiased branching ratio measurements, we also determine 
the branching ratio to the r~— 7r^7r^K+(7r'^) final state to ensure that the result is 
consistent with zero. In addition, we examine the resonant structure of these decays, 
and, under the assumption that the r~^z/^ K~7r~7r"'" final state is dominated by the Ki 
intermediate resonances, we determine the branching ratio of r~— Ki (1270) relative 
to the sum of the r~— Ki (1270) and r^— s^i/,- Ki(1400) branching ratios. 

The z/,- K~7r~7r"'"(7r'') final states occurring through r~— z/,- K~Kg(7r°) are con- 
sidered background in the analysis presented here. In addition, three-prong tau decays 
with charged kaons and more than one 7r° are severely phase-space suppressed and are 
therefore neglected. 



2 The OPAL Detector 

The OPAL detector is described in detail in reference The component of the 
detector most important to this analysis is the central jet chamber, which measures the 
momentum and specific ionisation energy loss, di^/dx, of charged particles. The jet 
chamber provides the only means of distinguishing between charged pions and kaons 
in the momentum range of interest to this study. 



4 



The jet chamber is a cyhnder 4 m long and 3.7 m in diameter, and is divided by 
cathode wire planes into 24 azimuthal sectors. Each sector contains one radial plane 
of anode wires, which are staggered to resolve left-right ambiguities. The chamber is 
contained in a solenoidal magnetic field of 0.435 T, and is filled with an argon-methane- 
isobutane gas mixture at a pressure of 4 atmospheres. This arrangement provides a 
transverse spatial resolution of cTxy = 130 /xm, and a two-hit resolution of 2.5 mm. In 
the barrel region of the jet chamber (| cos0| < 0.72) Q, the ionisation energy loss of a 
charged particle is sampled up to 159 times. A truncated mean is formed by discarding 
the highest 30% of the measurements, resulting in a d-E/dx resolution of about 3% for 
isolated tracks in the chamber 0. 

A layer of wire chambers with drift direction parallel to the z-axis lies immediately 
outside the jet chamber. These 2;-chambers accurately determine the polar angle of 
charged particles traversing the central detector. A lead-glass electromagnetic calorime- 
ter and presampler chambers are located outside the magnetic coil and jet chamber 
pressure vessel. The electromagnetic calorimeter is primarily used in this analysis to 
identify tau decays which include a vr" in the final state. The return yoke of the OPAL 
magnet is instrumented for hadron calorimetery, and is surrounded by external muon 
chambers. The hadron calorimeter and muon chambers are mainly used in this analy- 
sis to select control samples containing muons and to veto non-tau backgrounds in the 
tau-pair candidate sample. 



2.1 Parameterisation of Ionisation Energy Loss 

Figure |l|(a) shows the dependence of dE/dx on the momentum of various particles 
in the OPAL jet chamber, as predicted by the OPAL energy loss parameterisation, 
-Dprcd 0- Figure |I|(b) shows the particle resolving power, IZij, versus momentum for 
various pairs of particle species i and j, where 

Ti-ij = |-Dpj,g^ — -Dpi-cdl / ^1 

and where a is the mean of the uncertainties of -Dpred ^md -Dpj.gj. Note that the OPAL 
jet chamber yields a pion/kaon separation of at least 2cr for particles between about 2 
to 30 GeV/c. 

The parameters used in the calculation of Z^pred and a are tuned to the measured en- 
ergy loss of charged particles in hadronic Z° decays, yielding a d-E/dx parameterisation 
that is accurate enough for nearly all analyses of OPAL data. 

"'^A spherical coordinate system is used, with the +z-axis in the direction of the circulating electron 
beam. The angle is defined as the polar angle with respect to the +z-axis, and is defined as the 
azimuthal angle measured from the +a;-axis, which points towards the centre of the LEP ring. 
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However, hadronic tau decays are characterised by lower multiplicity and narrower 
jets than those produced in hadronic TP decays, resulting in small systematic differences 
between the measured dE/dx in the two environments. Although these differences are 
less than one a in magnitude, they must be corrected for accurate particle identification 
in three-prong tau decays. Based on a study of the measured energy loss of muons in 
e"'"e~— s> /i^/i~,r^— >/i^z7^ i/^, and e"'"e"— s>e"'"e~/i^/i~ samples, the following corrections are 
applied to the measured energy loss, -Dmeas, and to the OPAL dE jdx parameterisation 
and its uncertainty, Z^pred and a: 

Multiplicative correction to -Dprcd • To obtain a parameterisation which correctly pre- 
dicts the measured energy loss of particles from tau decays, D^j-^d must be cor- 
rected with a multiplicative factor, s{[3). Parameterising s{(3) as a second order 
polynomial in Q{[3) = — log(l — /?^) yields satisfactory results. The magnitude 
of the correction is of order 1% for tracks in one-prong tau decays. 

Additive correction to -Dmeas : The measured energy loss of charged particles in low- 
multiplicity events is found to depend strongly on the azimuthal separation, 0, 
between the track and anode plane in the jet chamber cell. To improve the 
dE/dx resolution, this behaviour is corrected with a function, /(</>), with seven 
parameters tuned to the measured energy loss of e"'"e~— *• candidates. The 

magnitude of this correction is also of order 1%. 

Multiplicative correction to a: To correctly predict the spread in the dE/dx of charged 
particles from one-prong tau decays, the OPAL parameterisation of a must be 
corrected with a multiplicative factor, Sres ~ 0.9. 

Although the application of these corrections significantly improves the parame- 
terisation of dE/dx in a low multiplicity environment, the energy loss distribution 
contains a significant non-Gaussian component. This feature is most apparent in the 
corrected dE/dx stretch distribution of the unlike-sign tracks of tau decays classified 
as r~— s-i/^ 7r~7r~X+(7r°) 0. Stretch energy loss. Si, is defined as: 

Si = (-Dmeas ^ -^pred)/ ' 

where i is the particle hypothesis used to calculate the predicted dE/dx and its uncer- 
tainty. The r~— >z/t- 7r~7r~X"'"(7r'') control sample is selected using Sk > 2.0 to identify 
three-prong tau events with two like-sign pion candidates in the final state. Since the 
r~— >z/t- 7r^7r^K+(7r°) final states are highly suppressed, the unlike-sign tracks in this 
sample are expected to be over 99.7% pure in pions. The stretch dE/dx distribution 
under a pion hypothesis of this sample is seen in Figure ||(a). 

^The unlike-sign track is the track in the three-prong decay with charge opposite that of the initial 
tau. 
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The non-Gaussian component is large in the region of the dE/dx distribution ex- 
pected to be populated by kaons. Parameterisation of the non-Gaussian component 
would therefore lead to very large statistical and systematic uncertainties in the mea- 
surement of the T~ ^ Ut K~ 7r~ 7r~^ (n^) and z/t- K~7r~K"'"(7r°) branching ratios. The 
tail is dramatically reduced in the dE/dx distribution of tracks in three-prong tau 
decays which lie closest to the anode plane, as shown in Figure 0(b). When tracks 
are spatially close together in the jet chamber cell, the tail of the pulse from the 
track closest to the anode will effectively shift the pedestal of the pulse of the second 
track, leading to an incorrect determination of deposited charge for that track, and 
the observed tails in the dE/dx distribution. The low-end non-Gaussian tail exists 
in the measured energy loss of charged particles in hadronic events, but is more 
pronounced in three-prong tau decays due to the tighter collimation of the particles in 
the decay. 

The simplest method to minimize this effect is to only consider the dE/dx of the 
track in the three-prong tau decay closest to the anode plane. Tracks meeting this 
requirement will be referred to as the 'first-track' three-prong sample. Further, the 
effect of pulses following the first is reduced by only considering first pulses separated 
from the second by more than 1 cm 0. 

The corrected stretch dE/dx distribution of unlike-sign tracks in the first-track 
r~— 7r~7r~X"'" sample, seen in Figure 0(b), is consistent with a Gaussian of zero 
mean and unit width. 



3 Monte Carlo Generated Event Samples 



For this analysis, tau lepton decays are simulated with the KORALZ 4.0 Monte 
Carlo generator and the Tauola 2.4 decay package To estimate biases in- 



troduced to the branching ratios by the tau-pair pre-selection procedure, we use a 
Monte Carlo sample of 600000 tau-pair events that have input branching ratios de- 
termined from world averages or theoretical expectations. In order to determine the 
branching ratios and to study the intermediate resonant structure of r~— z/^ K~7r~7r+ 
and r~— >z/^ K~7r~K'^ decays, large Monte Carlo samples of signal events are also pro- 
duced. All Monte Carlo samples are passed through a detailed simulation of the OPAL 



detector [|l^] and are subjected to the same analysis chain as the data. 



The T — >z/t-K tt 7r+ decay mode is generated in Tauola 2.4 primarily through the 
Ki(1400) intermediate resonance: 

T- z/^Ki(1400) ^ z/^K*(892)7r- ^ u^K-n-n+ (1) 

T" z/^Ki(1400) ^ z/^p(770)K" ^ z/^K-7r-7r+. (2) 
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However, the resonant structure of r — z/,- K vr 7r+ decays is likely more diverse than 
that represented in ([^) and (|), as the decay can also occur through the Ki(1270) inter- 
mediate resonance The Ki(1270) predominantly decays to K*(892)7r, KS(1430)7r, 
and p(770)K [|14|. To assess the sensitivity of the efficiency estimation to the assump- 
tions made about the r~— Ur K"?!"?!"*" intermediate resonant structure, a modified 
version of Tauola is used to generate the following decays: 

T- i/^Ki(1270) ^ z/^K*(892)7r" z/^K"7r"7r+ 

T- i/^Ki(1270) UrK*{U30)7i- UrK' 11-71+ 
T- z/^Ki(1270) z/^p(770)K- ^ z/^K-7r-7r+. 

In each case, the Ki(1270) decays through a Breit-Wigner resonance with an s- 
dependent width, and phase space limitations are taken into account |T3. Other 



Ki(1270) decay modes are ignored and interference effects between the various de- 
cay chains are neglected. Based upon experimental measurements of the Ki(1270) 
partial widths |T^, the Ki(1270) is assumed to decay 49% of the time to Kp, 32% of 
the time to K*(1430)7r, and 19% of the time to K*(892)7r. 

A sample of r^^z/^ K^tt^tt^ events is also generated through the Ki(1400) res- 
onance, again with the Ki decaying through a Breit-Wigner resonance with an s- 
dependent width. We only allow the Ki(1400)'s in this sample to decay to K*(892)7r, 
as experimental measurements indicate that this decay mode occurs much more fre- 



quently than the p(770)K mode [jT6|. 



The T — i>z/^K vr K"^ decay mode is generated in Tauola 2.4 through a mixture of 
the ai(1260) and p(1700) intermediate resonances: 

T- z/^ai(1260) z/^K*(892)K- z/^K"7r"K+ 
T- z/^p(1700) ^ z/^K*(892)K- u^K-n-K^^. 



Since the intermediate structure of the z/T-K~7r^7r+7r° and z/T-K~7r~K+7r'^ tau decay 
final states is not well understood theoretically, and experimental data on these de- 
cays are sparse, the Monte Carlo samples for these channels are generated through 
non-resonant production only. As a cross-check to ensure that the results of the 
r~— i>z/T- K~7r~7r+ and r~-^z/^ K~7r~K"'" resonant studies are independent of the mod- 
elling of the T--^i'rK~TT~TT'^TT^ and r~— >z/^ K~7r~K"'"7r° backgrounds, the following 
Monte Carlo samples are also used: 

T- z/^(K*(892)K7r)- ^ z/^K-tt-K+tt^. 



The second-order weak r"— > Ur 7r^7r^K+(7r°) decays are modelled using Monte Carlo 
samples of T~^i>r K~7r^7r+(7r°) decays, where the charge of the like-sign kaon and the 
unlike-sign pion are reversed. 
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The efficiency corrections in the branching ratio analysis are estimated using 
r^^z^T- K~7r~7r'^ and z/,- K~7r~K"*" events generated with the default version of 
Tauola 2.4, and r~— > K~7r^7r+7r° and T~—>-i'rK~7i~K~^7i^ events generated through 
non-resonant production. As a cross-check, the analysis is also performed with effi- 
ciency corrections derived from the r~ — >■ z/t-Ki(1270) and r~ — > z/,-Ki(1400) Monte 
Carlo samples. 



4 Branching Ratio Measurements 



4.1 Event Selection 

This analysis uses full data set collected by the OPAL detector between the years 
1990 and 1995 at e"'"e~ centre-of-mass energies close to the Z° resonance. Only data for 
which the tracking chambers and the electromagnetic calorimeter were fully operational 
are retained. The topology of e"'"e~— *-r"'"r~ events is characterised by a pair of back- 
to-back, narrow jets with low particle multiplicity. Jets are defined in this analysis 
by grouping tracks and electromagnetic clusters into cones with an opening half-angle 
of 35°, where each cone is assumed to contain the decay products of one of the tau 
leptons. Tau-pair candidates are selected by requiring two low-multiplicity jets with 
an average polar angle of |cos6'ave| < 0.90. Background events from other lepton 
pairs, hadronic Z" decays, and two-photon events are reduced with cuts on the event 
topology and total visible energy. These selections produce a sample of 147926 tau-pair 
candidates, with background = 2.73 ±0.27%, estimated by data control samples 

and Monte Carlo background samples. The tau-pair selection procedure is the same 



as that which is described in detail in reference [0 , except for the average polar angle 
selection, which is extended in this analysis to include tracks in the endcap region of 
the OPAL detector. 

The sum of the charges of the particles in each pre-selected tau-pair candidate event 
is required to be equal to zero. Three-prong tau decay candidates are selected from 
the jets in this sample containing three well-reconstructed charged tracks. The sum of 
the charges of the three tracks is required to be ±1, and decays with Kg candidates are 
excluded by rejecting tau decay candidates with one well-reconstructed neutral vertex 



with 7i~^'7i~ invariant mass between 0.4 and 0.6 GeV [0. The charge of the initial tau 
is taken to be the sum of the charges of the tracks in the cone associated with the tau 
lepton decay. This selection yields 38995 three-prong tau decay candidates. 
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4.1.1 Inclusive Candidate Samples (no 7r° Identification) 

The 'first-track' sample is formed by selecting the track in each three-prong tau decay 
candidate that is closest to the anode plane. This track is also required to lie within 
|cos^| < 0.90, and have at least 40 jet chamber hits used in the measurement of 
the dE/dx. All tracks in the first-track sample are also required to have reconstructed 
momentum between 3 and 90 GeV/ c. Tracks in the barrel region of the OPAL detector 
(I cos^^l < 0.72) are required to have at least 3 hits in the 2— chambers. Tracks which 
cross another track within the volume of the OPAL jet chamber are rejected from 
this sample. The 'like-sign' and 'unlike-sign' track samples are formed by selecting 
tracks in the first-track sample that have the same and opposite charge as that of the 
initial tau, respectively. The number of tracks in these two samples are 9296 and 4722, 
respectively. 

In first-order weak decays of the tau only the r^— z/,- K^7r^K+(7r'^) decay modes 
contribute kaons to the unlike-sign sample, whereas both the t^^Uj. K~7r^K+(7r°) and 
t^^Ut- K~7r~7r+(7r'') decay modes contribute kaons to the like-sign sample. Thus, after 
corrections for efficiencies and subtraction of background kaons, the numbers of kaons 
in these two samples are used to deduce the branching ratios to the Uj. K~7r"7r"*"(7r°) 
and i^T- K"7r~K"'"(7r°) final states. 

As a consistency check, we verify that the z/,- 7r~7r~K+(7r°) final states are 
indeed suppressed. To do this, we form a sample depleted in r^— > z/,- K~7r^K+(7r°) 
decays from the unlike-sign sample by requiring both associated like-sign tracks to 
have S'tt > —1.0. The K~ 7r~ -enhanced sample is formed from tracks in the unlike- 
sign sample which do no pass this selection. The number of tracks in these two samples 
is 1527 and 3195, respectively. The numbers of kaons in the i^'~7r~fr"'"-depleted sample, 
the i^^7r~ir"'"-enhanced sample, and the like-sign sample are used to simultaneously 
determine the tau branching fractions to the z/,- K~7r~7r"''(7r''), z/t- K~7r~K+(7r''), and 
Ur 7r~7r~K~^{n^) final states. 

4.1.2 7r° Identification 

To determine the branching ratios for decays with and without tt^'s, we divide all the 
inclusive samples into separate vr^-enhanced and vr^-depleted samples using J^E/ S-P > 
0.5 and J^E/ J2P ^ 0.5, respectively, where J^E/ J2P is the sum of the energy deposited 
in the electromagnetic calorimeter, divided by the scalar sum of the momentum of the 
charged particles in the tau decay cone. The J2E/J2P distributions for the Monte 
Carlo r~— >z/t- K~7r~7r+(7r°) and r~— i>z/^ K~7r~K"'"('7r°) samples are shown in Figure |^. 
The number of tracks in the vr'^-depleted samples with the like-sign, unlike-sign, and 
i^r^7r~iir"'"-depleted selections are 5984, and 3188, and 1110, respectively. 
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4.2 Estimation of the Number of Kaons in the Samples 



To determine the number of kaons in the data candidate samples described in 
Sections 4.1.1 and 4.1.2, we maximise a hkehhood function based on dE/dx: 

>C= n E h^''^ (3) 

j=l,N «=e,7r,K 

where W^^ is the dE/dx weight of charged particle j under particle hypothesis i, 

1 fDLs.-m)-m)Dt<i^ 



W'^ = . — exp 



(4) 



and where 



fi is the fraction of particle type i in the sample, where i is either pion, kaon, or 
electron. The values of fi are constrained to be non-negative, and the sum is 
constrained to be 1. 

N is the total number of particles in the sample, N = {Ne + A^tt + A^k)- 

-^meas mcasurcd dE/dx of the j^^ charged particle. 

-^pred is predicted dE/dx for charged particle j, calculated with the OPAL pa- 
rameterisation under particle hypothesis i , as derived from the measured dE/dx 
of charged particles in hadronic Z° events. 

(Tjj is the uncertainty on D^^^^^, calculated using the OPAL parameterisation, as de- 
rived from hadronic Z° events. 

Sres is the multiplicative correction to aij for all particle hypotheses i, and is allowed 
to vary without constraint in each fit. Separate Sres factors are used for charged 
particles in the barrel and endcap regions of the OPAL detector. The central 
values of these parameters are about 0.84 and 0.73, respectively. 

is the (3 dependent multiplicative correction to Dpred, and is equal to s{(3)+a, 
where s(/3) is determined from the one-prong control samples, and a is a cor- 
rection allowed to vary in each fit to compensate for possible slight differences 
between the dE/dx of the one-prong and three-prong decay environments. Sep- 
arate s'{(3) factors are needed for charged particles from three-prong tau decays 
in the barrel and endcap regions of the OPAL detector. These extra constant 
terms yield a correction to the dE/dx of approximately 1%. 

is the (f) dependent correction to the measured dE/dx, as obtained from the 
one-prong control samples. 
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Efforts are made to obtain a dE/dx parameterisation for tlie tau decay environment 
tliat is optimal for many particle species over a wide range of momenta. However, it is 
possible that the dE/dx corrections described in Section 2.1 may be somewhat more (or 
less) optimal for pions than they are for kaons in the momentum range of interest. Thus, 
to correct for any possible species-dependent quality differences in the parameterisation 
of dE/dx, an extra factor, Ck, is allowed to multiply the kaon predicted energy loss, 
D^j.^^, in the likelihood function found in Equation ^. This factor is determined in a 
likelihood fit to the measured dE/dx of one-prong tau decays to hadron final states to 
be Ck = 0.9968 ± 0.0014. 

Independent likelihood fits are performed for different ranges of momentum (13 bins 
of variable size between 3 and 90 GeV/c). A test of the fit with large Monte Carlo 
samples verifies that the resulting estimates for the kaon fraction have biases within the 
range — O.Scx/j^ to 0.3a at the 95% confidence level, where af^ is the typical statistical 
uncertainty on the kaon fraction returned by the fit. Thus, biases resulting from the 
fit procedure are neglected. 

The number of kaons found within each of the samples is given in Table ^, and 
the stretch dE/dx distributions of tracks in all momentum bins of the inclusive data 
like-sign and unlike-sign samples are shown in Figure |. The normalisation of the 
predicted distributions of the kaons, pions, and electrons in each sample is obtained 
from the results of the likelihood fits. The P^r degree of freedom between the first 
11 bins of the data and predicted distributions, which is the portion of the dE/dx 
distribution populated by kaons, is 15.9/10 and 10.6/7 for the like-sign and unlike-sign 
samples, respectively. 

4.3 dE/dx Systematic Studies 

A significant source of systematic uncertainty in this measurement is the uncertainty 
in the parameterisation of the predicted energy loss. To assess this systematic, we 
determine the sensitivity of the likelihood estimates of the number of kaons within 
each of the data samples to the uncertainties in the dE/dx correction factors obtained 
from the one-prong samples. To achieve this, we modify the likelihood function from 
equation ^ such that: 

where 

n is a vector containing the difference from the central values of the seven dE/dx 
correction factors that describe the function /(0), and the three correction factors 
that describe the function s(/5), where the central values are as derived from the 
fits to the one-prong control samples. 
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V is the covariance matrix for the dE/dx correction factors. 

In the first iteration, the correction factors are allowed to vary in the fit, and 
the returned values are found to be consistent with the input values. In the second 
iteration, the likelihood fit is repeated, this time keeping the dE/dx correction factors 
fixed to the values from the first iteration. The systematic uncertainty in /k due to 
the parameterisation of dE/dx is then obtained from the square root of the quadrature 
difference of the fit uncertainties in /k from the two iterations, and is shown in Table |^. 

4.4 Efficiency Correction 

The efficiencies for kaons from signal events in the pre-selected tau-pair sample 
to enter the candidate samples are estimated using signal events generated with the 
KORALZ 4.0 Monte Carlo generator and the Tauola 2.4 decay package, as described 
in Section 3. 

Table H shows the average efficiency estimates for the various signal channels in 
the pre-selected tau decay sample that contribute kaons to the data samples described 
in Sections 4.1.1 and 4.1.2. All efficiencies are corrected for biases introduced by the 
tau-pair pre-selection using a Monte Carlo tau-pair sample, and the efficiency uncer- 
tainties include the systematic uncertainty arising from this correction. The prese- 
lection biases are 0.935 ± 0.012, 0.925 ± 0.013, 0.908 ± 0.012, and 0.898 ± 0.013 for 
r^^z/^ K~7r^7r+,r^— > K^7r^K+,r^— K^7r^7r+7r°, and r^— K^7r~K+7r° decays, 
respectively. 

Figure ^ shows the efficiency versus momentum for z/^ K^7r^7r+ and 

T^^i'rK~iT~K~^ in the pre-selected tau-pair sample to contribute kaons to the in- 
clusive like-sign and unlike-sign samples, respectively. The branching ratio analysis is 
performed using an efficiency correction that is binned in momentum. 

The variation of the branching ratios due to alternative intermediate resonant 
structure scenarios for the r~— Ur K^tt^tt"'' final state is assessed using efficiency es- 
timates from signal events generated by the modified version of Tauola 2.4 through 
the Ki(1270) and Ki(1400) intermediate states as described in Section 3. The mix- 
ture of Ki(1270) and Ki(1400) is taken from the results of the analysis of the inter- 
mediate resonant structure of t~ K~ 7r~ candidates, as described in Section 5, 
R = Bi{t~^ z/^Ki(1270))/Br(r- ^ z/^ (Ki(1400) or Ki(1270))) = 0.71 ± 0.19. The 
branching ratios obtained within this range are found to be in agreement with the 
central values to within 1.5cr of the combined Monte Carlo statistical uncertainty. 

As a further cross-check, the branching ratios are evaluated using z/^ K^7r^7r+ 
and T^—>- Vr K~7r~K"'" efficiencies estimated using Monte Carlo events generated through 
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non-resonant production only. The resulting branching ratios are in agreement with 
the branching ratio central values to within la of the combined Monte Carlo statistical 
uncertainty. 

4.5 Kaon Background Correction 

Background kaons in the data samples described in Sections 4.1.1 and 4.1.2 are 
estimated by applying the same selection criteria to Monte Carlo samples of hadronic 
decays and tau-pair decays. From the number of selected events that contain kaons, 
the estimated number of background kaons are derived, as summarised in Table |^. 
Dominant sources of background kaons include low-multiplicity hadronic Z° events, 
and — > Ut-K" > neutrals decays. 

4.6 The Branching Ratio Calculation 

The branching ratios in the exclusive tau decay channels of interest are calculated 
from the numbers of kaons, as estimated by the likelihood fit, in several momentum 
bins of the exclusive candidate samples, listed in Table ^ The number of kaons within 
each momentum bin i of candidate sample j is corrected for background as described 
in Section 4.5, yielding 

where 

jypre-sei jg numbcr of pre-selected tau decay candidates. There were 295852 tau 
decay candidates recorded in the OPAL detector between the years 1990 and 
1995. 

ynon-T -g ^^le estimated background from non-r sources in the pre-selected tau decay 
candidates {f^°^-^ = 2.73 ± 0.27%). 

A^^ is the number of kaons in momentum bin i of candidate sample j, as estimated 
by the likelihood fit. The number of kaons summed over all momentum bins of 
each sample is shown in Table 0. 

^bkgnd is estimated number of background kaons in momentum bin i of candidate 
sample j, as estimated from Monte Carlo tau decay and hadronic Z" decay sam- 
ples. The number of background kaons summed over all momentum bins of each 
sample is shown in Table 0. 
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A set of linear equations that relate the Rij to the tau branching ratios are solved 
simultaneously to determine the branching ratios. The set of equations includes effi- 
ciency corrections as a function of momentum for each decay channel, as determined 
from Monte Carlo generated events, as described in Section 4.4. The efficiency correc- 
tion does not assume that the momentum spectrum of the signal kaons in each sample 
follows those of kaons from the Monte Carlo. For background corrections, however, the 
momentum spectra of Monte Carlo generated events are used. 

In the first step of the branching ratio calculation, the r"— i>z/^ 7r~7r~K^(7r'^) branch- 
ing ratios are not assumed to be zero. The branching ratio calculation uses the division 
of the unlike-sign samples into the K~'7r~K"'"-enhanced and K~7r~K"'"-depleted samples 
in conjunction with the like-sign samples, for a total of six exclusive samples, as listed 
in Table ^ The results of the calculation are: 
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where the uncertainties are statistical only, and where the inclusive branching ratios 
are obtained from the sums of the exclusive branching ratios. Correlations between 
the exclusive branching ratios are taken into account in this calculation. The results 
for the branching ratios to the z/,- 7r~7r^K+(7r'^) final states, which can only occur as 
second-order weak interactions, are consistent with zero, as expected, indicating that 
there are no significant biases in the kaon identification procedure. 

To obtain more precise branching ratios for the first-order weak decays, we set 
the r~— > z/t- 7r^7r^K+(7r'^) branching ratios to zero and repeat the calculation. In this 
iteration, the calculation uses the unlike-sign and like-sign exclusive samples, for a total 
of four samples. The results of this calculation are: 

Br(r- ^ z/^ K-7r-7r+(7r°)) = +0.343 ± 0.073 ± 0.031 % 
Br(r" z/^K"7r"K+(7r°)) = +0.159 ib 0.053 =b 0.020 % 

Br(r- z/^ K-7r-7r+7r°) = -0.017 + 0.076 + 0.060 % 
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Bi{t- Ur K-7r-K+7r°) = +0.072 ± 0.085 ± 0.051 % 

Br(r- ^ i^r K-n-n+) = +0.360 ± 0.082 ± 0.048 % 
Br(r- ^ z/^ K-7r-K+) = +0.087 + 0.056 + 0.040 %, 

where the first uncertainties are statistical and the second are systematic. The sum- 
mary of the systematic uncertainties for each branching ratio is shown in Table ^. The 
correlation matrix for the exclusive branching ratios is found in Table 

The z/t- K^7r^7r+(7r°) branching ratios from the first and second calculations 
are approximately 70% correlated, and the t^^v^ K^7r^K+(7r°) branching ratios from 
the first and second calculations are approximately 50% correlated. In both cases, the 
resulting differences in the branching ratios between the two iterations are consistent 
to within 1.3cr of the combined statistical uncertainties. 

The central values of the branching ratios are evaluated using an ^E/ Y^P selection 
of 0.5 to distinguish between states including and not including a 7r°. The branching 
ratios are also evaluated using J^E/ J2P selections of 0.3, 0.4, 0.6, and 0.7. The branch- 
ing ratio systematic uncertainty associated with the J2P selection is taken as the 
RMS spread of the five values, and is quoted in Table ^. 



5 Resonant Structure 



5.1 Event Selection 



The selection of events used in the study of the resonant structure begins with the 
three-prong tau decay candidate sample described in Section 4. Events which include 
an e"'"e~ pair from a photon converting in the detector material are identified and 
rejected on the basis of a topological conversion finder |T^. All three tracks in each 



decay are required to have at least 40 jet chamber hits used in the measurement of 
dE/dx. Tracks in the barrel region of the OPAL detector are required to have at least 
3 hits in the 2;-chambers. 

Using 5^ < —2.0 to identify kaons and 15*7^1 < 1.5 to identify pions, events are 
classified into the following channels: 

• t^^Ut- K^tt^tt^ candidates have one like-sign kaon, and an unlike-sign and like- 
sign pion. 

• T^^Ur K^TT^K^ candidates have one like-sign pion, and an unlike-sign and like- 
sign kaon. 
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• T 71 7r+(7r°) candidates have all three tracks identified as pions. 

The number of kaons in the lowest momentum bin of the candidate samples used in 
the branching ratio analysis is found to be consistent with zero. Thus, to increase the 
signal to noise ratio in the samples used in the resonant studies, all kaon candidates are 
required to have momentum of at least 5 GeV/c. Pion candidates are required to have 
momentum of at least 0.75 GeV/c. The above selections produce 950 t^^v^ K^7r^7r+ 
and 79 r^^z/^ K~7r~K"*" candidates. 

The non-Gaussian tails in the data dE /dx distribution are not as problematic in this 
portion of the analysis as they are in the branching ratio analysis and thus no attempt is 
made to remove them. The differing shapes of the signal and the t^^z^t- 7r^7r^7r+(7r'^) 
invariant mass distributions yield some separation power between signal and back- 
ground in the data samples, and the dE /dx selections used to identify pions and kaons 
in this study are designed only to provide three-prong tau decay samples enhanced 
with r^— > z/^ K^7r^7r+ or r^^z/,- K^tt^K"*" decays. Other than these initial selections, 
dE /dx is not used to distinguish between signal and background. 



5.2 r —^Ut-K tt 7r+ 

The background in the r^^z/,- K^7r^7r+ sample consists primarily of 
t^—^Ut- 71^ 71^71^(71^) events, along with contributions from r~— z/,- K^tt^tt+tt'^, other 
tau decays, and hadronic Z° events. 

To analyse the resonant structure of these decays, we examine the Mi^-^-^+, 
MK-7r+, and Mtj-tj+ invariant mass distributions of the candidates. To estimate the 
shape of the dominant three-pion background in these distributions, we use the data 
r~— >z/t- 7r"7r^7r+(7r°) candidate sample, with the 5 GeV/c momentum selection placed 
on the particle which corresponds to the like-sign kaon candidate. This selection pro- 
duces a sample of about 16000 events. The shape of the r~— z/^ K^Tr^vr+vr^ background 
in the distributions is parameterised with Monte Carlo r^^z/,- K^tt^tt+tt^ events gen- 
erated through non-resonant production. Other background fractions are expected to 
be on the order of a percent or less and are neglected. 

The description of the shape of the signal portion of these invariant mass distri- 
butions depends upon assumptions made about the intermediate resonant structure of 
T^^UrK^n^n'^ decays. We consider four different scenarios: 

1. The T^^Uj. K^7r^7r+ decays occur through a mixture of Ki(1270) and Ki(1400). 
For this study we assume that the Ki(1270) has a width of 90 MeV, and the 
Ki(1400) has a width of 174 MeV, the current world average widths of these 
resonances IH . 
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2. The r^— »z/^ K^7r^7r+ decays occur through the Ki(1270) and Ki(1400) reso- 
nances, and that the widths of these resonances are both 300 MeV, as suggested 
in reference 0. 

3. The r^^z/,- K^7r^7r+ decays occur through the K*(892)7r~ and p{770)K~ inter- 
mediate resonances. 

4. The r"— K~7r~7r+ decays occur through non-resonant production (phase 
space) only. 

Monte Carlo samples having at least 10 times the statistics of the expected signal com- 
ponent of the data T^^UrK^n^n^ candidates are generated with a modified version 
of Tauola 2.4 under each of these resonant structure assumptions. 

To estimate the background and signal fractions in a way which accounts for corre- 
lations between the invariant mass distributions, we divide the MK-7r-7r+, MK-7r+, and 
M^-^+ invariant masses of the data, signal Monte Carlo sets, and background samples 
into 7 bins each, to form a 7x7x7 matrix. Further discrimination between signal 
and background is obtained by dividing the scalar sum of the momenta in each three- 
prong decay into 24 bins between and 48 GeV/c. The correlations between J2 P and 
the invariant masses are not significant compared to the inter-correlations between the 
invariant masses. 

Using the various signal and background distributions as templates, a binned maxi- 
mum likelihood fit is performed simultaneously to the data 7x7 x7 matrix and the P 
distribution to determine the most probable fractions of signal and r~— z/^ K^tt^tt+vt'' 
background in the r~— K~'7r~7r+ candidates. For the first three scenarios, two signal 
fractions are allowed to float. Bins in the matrix and Y1,P distribution which contain 
data events but no predicted signal or background events are neglected in the fit 0. A 
test of the procedure using a Monte Carlo event sample of the same size and approxi- 
mate composition of the data sample reveals no significant bias in the estimates of the 
signal and background fractions from the fit. 

The results of the fits are shown in Table ^. Table ^ also gives the per degree of 
freedom between the data and predicted MK-,r-7r+, ^K-7r+, and M^-^+ distribu- 
tions, where the normalisation of the predicted distributions is obtained from the fit 
results. The binning of these distributions is the same as that seen in Figure ^ which 
displays the results of the fit to the data using the assumption of the world average 
widths for the two Ki resonances. The correlation between the estimated fractions of 
Ki(1270) and Ki(1400) obtained from this fit is about —0.30. A fit where the signal 
and z/^ K~7r^7r+7r° fractions are set to zero yields — 21og£/£max = 107. Since this 

•^Only one such bin exists in the fit for the central values of the signal and background fractions, 
and it contains only one event. 
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fit has three fewer degrees of freedom than the first, the difference in — 2Iog£ should 
be distributed with three degrees of freedom if the sample does indeed consist of 
only three-pion background. This test disfavours this hypothesis at a confidence level 
of over 99%. 

Table |^ includes estimates of the efficiencies for z/t- K^tt^tt^ and 

T~^h'rK~TT'7r~^TT^ in the pre-selected tau-pair sample to contribute to the 
r^— s> z/t- K^TT^vr^ candidate samples. These efficiencies are not corrected for biases 
introduced by the r~— > z/,- K'tt"??"*" selection procedure. In order to verify that the 
signal fractions returned by the fit are reasonable, we use these efficiencies and the fit 
fractions to calculate estimates of the r^^z/^ K^7r^7r+ and r^— s>z/^ K^tt^tt+tt^ branch- 
ing ratios, which are also given in Table |^. Although the efficiencies have not been 
corrected for bias, these estimates are in agreement with the branching ratios obtained 
in Section 4 for all resonant structure assumptions, except for r~— ^z/,- K~n~'ir~^ through 
phase space only. 

As a cross-check, Monte Carlo generated — >• z/t-K^cc; decays are used to es- 
timate the shape of the T~—>-i'rK~'ir~7i'~^Tr^ in the data distributions, rather than 
r^-^ z/t- K^tt^tt+tt'^ generated through non-resonant production. For all four assump- 
tions made about the r~— > z/,- K^tt^tt^ intermediate resonant structure, these fits re- 
turn signal fractions within 0.15(7 of the central values, and r~^z/^ K^tt^tt+tt^ fractions 
within 0.35(7 of the central values, where a refers to the statistical uncertainties of the 
central values of the fractions. 



5.2.1 r- z/^Ki 

Under the assumption that resonant structure of t~^Vj. K~7r^7r+ decays is dominated 
by the Ki intermediate resonances, we derive 

^ fraction of Ki(1270) in K^tt^tt"^ candidates 69 ± 16 

~ fraction of Ki(1400)+Ki(1270) in K-ji-ti^ candidates 

where the uncertainty is statistical, and where the fractions of Ki(1270) and Ki(1400) 
are taken from the fit which assumes the widths of these resonances are 90 MeV and 
174 MeV, respectively. 

The signal and background template distributions used in the binned likelihood 
fit have statistical uncertainties associated with them. To determine the systematic 
uncertainties on the fit fractions arising from these uncertainties, the fit is repeated 25 
times, each time randomly varying each bin of the template distributions by sampling 
Poisson distributions with means equal to the original bin contents. The resulting 
systematic uncertainties are included in Table ||, and are derived from the RMS spread 
of the signal fraction estimates from the 25 fits. 
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To check for undue variation that may be produced by the choice of binning used 

in the fit, the above procedure is repeated for various binning schemes. In all cases, 
the RMS variation of the results returned by the different fits is less than the statis- 
tical uncertainty from the original fit. The variation of i?fit due to uncertainty in the 
Ki(1270) branching fractions is also studied, and found to be negligible. 

The value of i?fit is corrected for the efficiencies for Ki(1270) and Ki(1400) in the 
pre-selected tau-pair sample to contribute to the Vj. K"?!"?!"*" candidate sample 
(0.052 ± 0.001 and 0.058 ± 0.002, respectively). This yields: 

Br(r-^ ..K.(1270)) ^ 0.71 ± 0.16 ± 0.11, 

Br(T- ^ z/,Ki(1400)) + Br(r- ^ z/,Ki(1270)) 

where the first uncertainty is due to the statistical uncertainty from the fit, and the 
second arises from the limited statistics of the Monte Carlo generated samples and the 
contribution from other systematic effects. 

To determine if R depends on assumptions made about the widths of the Ki reso- 
nances, the above procedure is repeated, simulating the signal portion of the distribu- 
tions using r~— > K~7r~7r"'" Monte Carlo events generated through the Ki(1270) and 
Ki(1400) resonances, both with width 300 MeV. This procedure yields: 

R = 0.68 ±0.13 ±0.11, 

which is in agreement with the central value of R to within 0.2(7 of the combined Monte 
Carlo statistical uncertainty. 



5.3 r-^zy^K-7r-K+ 

The T~— K~7r~K+ candidate sample is expected to consist of t~^VtT^~t^''t^^{'k^) 
background, along with r^^z/^ K~7r~K+ signal, some t~— ^z/,- K^7r^K+7r° contamina- 
tion, and some contamination from other tau decays and hadronic Z° events. 

To analyse the resonant structure of t~— ^i^^ K^tt^K^ decays, we examine the 
MK-7r-K+, -^K-K+, and M7r-K+ invariant mass distributions of the candidates. To 
estimate the shape of the three-pion background in these distributions, we use the 
data t~^Vt-tt~tt~tt^{t[^) candidate sample, with the 5 GeV/c momentum selection 
placed on the particles which correspond to the two kaon candidates. This selection 
produces a sample of over 10000 events. 

The shape of the r~— > z/.^ K~7r^K+7r'^ background in the distributions is parame- 
terised with Monte Carlo t~— ^i/^ K~7r~K"'"7r° events generated through phase space. 
Other background fractions are expected to be on the order of a percent or less and 
are neglected. 
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To describe the shape of the signal portion of these invariant mass distributions, 
we consider three different scenarios: 

1. The r~— > z/^ K~7r~K"'" decays occur through the same intermediate resonant struc- 
ture used in the default version of Tauola 2.4, as described in Section 3. 

2. The T'^u^K-n-K+ decays occur through K*(892)K-. 

3. The T^^UrK^TT^K^ decays occur through phase space only. 

Monte Carlo samples having at least 10 times the statistics of the expected signal 
component of the data t^^Ut-K^ti^h^ candidates are generated under each of these 
resonant structure assumptions. 

To determine the composition of the candidate sample, we divide the MK-7r-K+; 
Mk-k+; and Mjr-K+ invariant masses of the data, signal Monte Carlo, and background 
samples into 24 bins each. Further discrimination between signal and background is 
obtained by dividing the scalar sum of the momenta in each three-prong decay into 24 
bins between and 48 GeV/c. Since this study is qualitative in nature only, correlations 
between the distributions are ignored. 

The results of the fits using the three resonant structure assumptions are shown 
in Table ||. Table ^ also gives the per degree of freedom between the data and 
predicted J^P, ^K-7r-K+) ^k-k+; and M^-k+ distributions, where the normalisation 
of the predicted distributions is obtained from the fit results. The binning of these 
distributions is the same as that seen in Figure |^, which displays the results of the best 
fit to the data. A fit where the signal and z/^ K^tt^K+tt^ fractions are set to zero 
yields — 21og£/£max = 51. Since this fit has two fewer degrees of freedom than the 
first, the difference in — 21og£ should be distributed with two degrees of freedom 
if the sample does indeed consist of only three-pion background. This test disfavours 
this hypothesis at a confidence level of over 99%. 

Table ^ includes estimates of the efficiencies for r^^z/^ K^7r^K+ and 
r^^z/^ K~7r~K"^7r° in the tau-pair sample to contribute to the r~— ^z/,- K~7r~K"*" can- 
didate samples. These efficiencies are not corrected for biases introduced by the 
z/t- K~7r^K+ selection procedure. In order to verify that the signal fractions re- 
turned by the fit are reasonable, we use these efficiencies and the fit fractions to calcu- 
late estimates of the r~— s>z/^ K~7r~K^ and r~— >z/,- K~7r~K^7r° branching ratios, which 
are also given in Table |^. Although the efficiencies have not been corrected for bias, 
these estimates are in agreement with the branching ratios obtained in Section 4 for 
all resonant structure assumptions. 

As a cross-check, Monte Carlo generated r~— >z/T-(K*(892)K7r)^ decays are used to 
estimate the shape of the r^^z/,- K^tt^K+tt'^ in the data distributions, rather than 
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r^— >z/t- K^7r^K+7r° generated through non-resonant production. For all three assump- 
tions made about the r~— >z/^ K~7r~K"^ intermediate resonant structure, these fits re- 
turn signal fractions within 0.25a of the central values, and r~— > z/,- K~7r^K+7r° frac- 
tions within 0.40cr of the central values, where a refers to the statistical uncertainties 
of the central values of the fractions. 



6 Summary and Discussion 



From a sample of 295852 tau decays recorded in the OPAL detector between the 
years 1990 and 1995, we determine the branching ratios: 



Br(r' i/^K-7r-7r+(7r°)) = 0.343 ± 0.073 ± 0.031 



'0 



•'0 



Br(r- ^ z/^K-7r-K+(7r°)) = 0.159 ± 0.053 ± 0.020 

Br(r- ^ z/, K-7r-7r+7r°) < 0.179% (95% CL) 
Br(r- z/,K-7r-K+7r°) < 0.247 % (95% CL) 

Br(r- ^ z/^ K-7r-7r+) = 0.360 ± 0.082 ± 0.048 % 
Bi{r- ^ UrK-n-K+) = 0.087 ± 0.056 ± 0.040 %, 

where the first uncertainty is statistical and the second is systematic. 

The r~— > Ur K~7i~7r~^ inclusive and exclusive branching ratios are in agreement with 
both theory and previous empirical measurements, as listed in Table |1]. 

In a separate analysis, we explore the resonant structure of t^—^Ut-K^ti^tt^ and 
z/^ K~7r^K+ decays. Under the assumption that the intermediate resonant struc- 
ture of the tau decay to the u^- K~7r~7r"'" final state is dominated by the Ki intermediate 
resonances, we determine: 

Br(r- ^ z/^Ki(1270)) 
Br(r- ^ z..K,(1400)) + Br(r- ^ z..Ki(1270)) " O-^l ± 0-16 ± O.IL 



There are two previously published results for R: 



R = 0.35Zlil TPC/27 1994 I 
R = 0.91 ± 0.29 ALEPH 1999 M 



with an average of 0.63 ± 0.21. The OPAL result is in agreement with this average. 
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It has been suggested in that theoretical predictions best match the world av- 
erages for the r~ —>■ (Ktttt)^ branching ratios if the Ki resonances are in fact wider 
than the current world average widths, rKi(i270) = 90 MeV and rKi(i4oo) = 

174 MeV. 

From SU{3) flavour symmetry arguments, reference suggests that the actual widths 
of these resonances are likely to be approximately the ai(1260) width (250 MeV or 
greater). We find that the data do indeed favour wider Ki resonances, but that the 
world average widths are also consistent with the data. 
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T- DECAY 
MODE 


Strangeness 


EXPERIMENT 
BR (%) 


THEORY 

BR (%) 


Vr K^7r^7r+(7r°) 


-1 


0.275 ±0.064 ALEPH98lT|] 

0.421 ±0.068 CLEO980] 
0.58i|]:ll TPC/2794§ 

0.22i[J:l^ DELC085|] 
0.343 ± 0.079 (this analysis) 




z/^K-7r-K+(7r°) 





0.238 ± 0.042 ALEPH98ll|] 

0.178 ±0.036 CLE098||] 
0.151°:°^ TPC/2794[| 
0.159 ± 0.057 (this analysis) 


- 




-1 


0.214 ±0.047 ALEPH98ll|] 
0.346 ±0.061 CLEO980] 
0.360 ± 0.095 (this analysis) 


0.77 i 

0.35 to 0.45 [|] 

0.18 e 


K-7r-K+ 





0.163 ±0.027 ALEPH98ll|] 

145 ±0 031 CLE098i21 
0.22i|]:}^ DELC085|] 
0.087 ± 0.069 (this analysis) 


0.20 i 
26 W\ 




-1 


0.061 ± 0.043 ALEPH98ll|] 
0.075 ±0.032 CLE098||] 
< 0.179 (this analysis) 




z/^ K-tt-K+ttO 





0.075 ±0.033 ALEPH981T|] 
0.033 ±0.019 CLE098||] 
< 0.247 (this analysis) 


- 


Vr 7r^7r^K+(7r°) 


+1 


< 0.25 TPC/2794[| 




z/^K-K-7r+(7r°) 


-2 


< 0.09 TPC/2794[| 




i/^K-K-K+(7r°) 


-1 


< 0.21 TPC/2794[| 

< 0.019 ALEPH981T|] 





Table 1: Three-prong decays of the tau lepton which include charged kaons in the 
final state (charge conjugate decays are implied). Experimental uncertainties are the 
combined statistical and systematic uncertainties. Branching ratios quoted as limits 
are the 95% confidence limits, and the (vr") notation refers to decay modes with or 
without an accompanying 7r°. 
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TT^-depleted Sample 


7r"-enhanced Sample 


like-sign 


unlike-sign 


like-sign 


unlike-sign 


K-7r-K+ 
enhanced 


K-7r-K+ 
depleted 


K-7r-K+ 
enhanced 


K-7r-K+ 
depleted 


iVK 
-^bkgnd 


128.4 ± 16.5 ±2.3 
5.8 ± 1.8 


36.7 ± 9.5 ± 1.2 



1.4 ±3.7 ±0.5 




29.5 ± 9.2 ± 1.1 
8.5 ±2.1 


5.9 ±5.6 ±0.6 
2.5 ± 1.5 


4.6 ±3.7 ±0.3 



fK-7r-K+ 
^■K-n-K+ 

^K~7r~7r+7rO 
CK-TT-K+TrO 
^7r~7r~K+7rO 


0.094 ± 0.003 
0.078 ± 0.003 

< 0.001 

0.044 ± 0.002 
0.048 ± 0.003 

< 0.001 


< 0.001 
0.104 ± 0.003 
0.065 ± 0.006 

< 0.001 
0.053 ±0.003 
0.033 ±0.002 


< 0.001 
0.003 ±0.001 
0.029 ± 0.002 

< 0.001 
0.001 ±0.001 
0.012 ±0.001 


0.013 ± 0.001 
0.007 ±0.001 

< 0.001 

0.047 ±0.002 
0.030 ±0.002 

< 0.001 


< 0.001 
0.008 ±0.001 
0.008 ±0.001 

< 0.001 
0.032 ±0.002 
0.036 ±0.001 


< 0.001 

< 0.001 
0.005 ± 0.001 

< 0.001 

< 0.001 
0.011 ±0.001 



Table 2: The number of kaons in each data sample, A^k, where the first uncertainty is the statistical uncertainty from the 
maximum likelihood fit to the measured dE/dx of the tracks in the sample, and the second is the systematic uncertainty 
arising from the imcertainties in the parameterisation of dE/dx. Also shown are the estimated backgrounds and the average 
efficiencies for the various signal channels which contribute to each sample. 





Branching Ratios (%) 


-DK-7r-7r+(7rO) 


^K-7r-K+(7rO) 


JJK-7r-7r+7rO 


JJK-7r-K+7rO 


JJK-7r-7r+ 


JJK-7r-K+ 


Central value 


0.343 


0.159 


-0.017 


0.072 


0.360 


0.087 


a (stat) 


±0.073 


±0.053 


±0.076 


±0.085 


±0.082 


±0.056 


a {dEjdx sys) 


±0.017 


±0.010 


±0.020 


±0.015 


±0.013 


±0.016 


a (MC Stat) 


±0.023 


±0.012 


±0.029 


±0.018 


±0.029 


±0.018 


(7(E^/E^sys) 


±0.012 


±0.012 


±0.049 


±0.045 


±0.036 


±0.032 



Table 3: Summary of the branching ratio central values and sources of uncertainty. 





BK-7r-7r+ 


BK-7r-K+ 


Bk-tt-tt+ttO 


BK-TT-K+Tr" 


BK-7r-7r+ 


+ 1.0 


-0.3 


-0.6 


+0.1 


BK-7r-K+ 


-0.3 


+ 1.0 


+0.4 


-0.8 


Bk-tt-tt+ttO 


-0.6 


+0.3 


+1.0 


-0.5 


Bk-tt-K+ttO 


+0.1 


-0.8 


-0.5 


+1.0 



Table 4: Correlations between the exclusive branching ratios. 
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Ki(1270) 

and Ki(1400) 


Ki(1270) 
and Ki(1400) 


K*(892)7r- 
and p(770)i^- 


Phase Space 




ri27o = 90 MeV 


ri27o = 300 MeV 








ri4oo = 174 MeV 


Tuoo = 300 MeV 






/'signal#l 


/ki(1400) = 0.043 ± 0.025 ± 0.017 


/ki(1400) = 0.065 ± 0.025 ± 0.021 


/k^tt = 0.100 ±0.027 


0.029 ± 0.033 


/signal#2 


/ki(i270) = 0.096 ± 0.034 ± 0.028 


/ki(1270) = 0.121 ± 0.039 ± 0.036 


fpK = 0.040 ± 0.235 






0.066 ± 0.028 


0.048 ± 0.036 


0.076 ± 0.030 


0.082 ± 0.029 


-2 log £/ £max 


8 





3 


96 


(EP) 

(Mk-.-.+ ) 
(Mk-.+ ) 


15.9/18 
26.5/20 
22.3/17 
26.1/18 


17.4/18 
23.7/20 
19.5/17 
27.8/18 


17.3/18 
26.7/20 
21.2/17 
26.3/18 


22.9/18 
29.0/20 
37.0/17 
25.9/18 




0.22 


0.22 


0.22 


0.21 


*^K~7r~7r+7rO 


0.17 


0.17 


0.17 


0.17 


BK-7r-7r+(%) 


0.207 ±0.054 


0.271 ± 0.058 


0.205 ± 0.334 


0.045 ± 0.052 


BK-7r-7r+7ro(%) 


0.128 ±0.054 


0.093 ± 0.069 


0.147 ±0.058 


0.158 ±0.056 



Table 5: Summary of the r~— > z/t- K~7r~7r"'" signal and r~— z/^- K~7r~7r"'"7r'^ background fractions obtained from the binned 
maximum likelihood fits to the data r~— > i^r K~7r~7r~^ candidate invariant mass and J2 P distributions. The first uncertainties 
are the statistical uncertainties from the fits. The second uncertainties for /ki(i4oo) and /ki(i27o) a-re the systematic 
uncertainties due to the limited statistics in the Monte Carlo generated samples. Also shown are the values of —2 log >C/>Cinax 
from the fit (relative to the best fit), and the resulting x^ per degree of freedom between the various data and predicted 
distributions. Estimates of the efficiencies for signal and r~— > z/^ K~7r~7r+7r° decays in the tau-pair candidate sample to 
contribute to the r~— > K~7r~7r"'' sample are given. From these efficiencies and the fraction estimates from the fits, 
estimates of the branching ratios are calculated. These branching ratios have not been corrected for biases introduced 
by the event selection and are used as a cross-check of the analysis only. The branching ratio uncertainties are from the 
statistical uncertainties of the fit fractions. 





Tauola 2.4 default 








mix of ai(1260) 
and p(1700) 


K*(892)K- 


Phase Space 


/K-7r-K+ 


0.473 ± 0.057 


0.447 ± 0.074 


0.517 ±0.086 


/k-tt-K+ttO 


0.000 ±0.304 


0.141 ±0.059 


0.109 ±0.061 


-2 log C/Craax 
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(Mk-^-k+) 

(Mk-k+) 

(M^-K+) 


7.6/5 
7.0/6 
5.5/6 
3.8/5 


7.3/5 
8.7/6 
10.6/6 
4.5/5 


8.4/5 
7.8/6 
8.5/6 
8.6/5 


CK-7r-K+ 


0.13 


0.14 


0.15 


%-7r-K+7rO 


0.10 


0.10 


0.10 


BK-7r-K+(%) 


0.096 ±0.012 


0.085 ±0.014 


0.096 ±0.016 


BK-7r-K+7ro(%) 


0.000 ± 0.079 


0.037 ±0.015 


0.028 ±0.016 



Table 6: Summary of the r~— > z/^ K~7r^K^ signal and r^^z/^ K^tt^K^tt back- 
ground fractions obtained from the binned maximum likelihood fits to the data 
t~^i/t- K~7r~K"'" candidate invariant mass and J2P distributions. The uncertain- 
ties are the statistical uncertainties from the fits. Also shown are the values of 
— 21og£/£max from the fit (relative to the best fit), and the resulting per degree 
of freedom between the various data and predicted distributions. Estimates of the 
efficiencies for signal and r~— > z/,- K^tt^K^tt^ decays in the tau-pair candidate sample 
to contribute to the T~^Vr K~7r~K"'' sample are given. From these efficiencies and the 
fraction estimates from the fits, estimates of the branching ratios are calculated. These 
branching ratios have not been corrected for biases introduced by the event selection 
and are used as a cross-check of the analysis only. The branching ratio uncertainties 
are from the statistical uncertainties of the fit fractions. 
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OPAL 




Momentum (GeV/c) 

Figure 1: (a) shows the ionisation energy loss i^pred ^ a function of the momentum 
for various particles in the OPAL jet chamber, (b) shows the resolution power TZy 
expressed in terms of the dE/dx resolution a, for various pairs of particle species i and 
J- 
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Figure 2: (a) is the stretch dE/dx distribution under the pion hypothesis for the 
unhke-sign tracks in the data z/^ 7r~7r^X+(7r*^) candidate sample (points). This 
sample is expected to have less than three kaons present. The curve is a unit width 
Gaussian of zero mean, whose normalisation is fit to the central part of the distribution. 
A significant non-Gaussian component is evident, (b) is the same distribution for the 
unlike-sign tracks in the data first-track t~— 7r~7r~X"'"(7r'^) candidate sample. 
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Figure 3: (a) shows the Z^-B/ Z^P distributions for Monte Carlo generated 
r~— > K^Tr^Tr"*" and t~— K~7r~7r"'"7r°. (b) shows the same distributions for Monte 
Carlo generated r^^i/^ K~7r~K+ and t~— ^i/,- K~7r~K+7r°. 
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OPAL 




Figure 4: (a) shows the stretch dE/dx distribution of tracks in the data inclusive 
hke-sign sample (points). The overlaid curves are the predicted distributions for the 
kaons, pions and electrons in the sample, assuming a Gaussian resolution function for 
the dE/dx measurements. The normalisation of the curves comes from the results 
of the likelihood fit. (b) shows the same distribution for tracks in the data inclusive 
unlike-sign sample. 
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Figure 5: (a) shows the efficiency as a function of kaon momentum for Uj. K~7r~7r+ 
in the pre-selected tau-pair sample to contribute kaons to the hke-sign sample; the solid 
points represent the default Tauola 2.4 generation of t~— ^i/^- K~7r~7r+ and the dotted 
points are calculated using r~— > i/t K^Tr^Tr"^ generated through phase space only, (b) 
shows the efficiency versus kaon momentum for t~— K~7r~K+ in the pre-selected 
tau-pair sample to contribute kaons to the unlike-sign sample; the solid points represent 
the default Tauola 2.4 generation of T~^i>r K~7r~K+ and the dotted points are calcu- 
lated using T'^Vr K~7r~K+ generated through phase space only. The branching ratio 
central values are calculated using K^tt^tt"'" and t~— ^i/,- K~7r~K'*' efficiencies 

estimated with the default version of Tauola 2.4. 
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OPAL t'^vK'tt'ti^ candidates 
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Figure 6: (a) is the J2P distribution and (b), (c), and (d) are the invariant mass 
distributions of the data T^—^UrK^n^Ti^ candidate sample (points). The histogram 
represents the predicted distribution, in which the normahsation of the background and 
signal components comes from the best fit results of the binned maximum likelihood 
fit which took the correlations between the invariant mass distributions into account. 
The shape of the signal is estimated using t~^i't-K~tt~tt^ events generated through 
the Ki(1270) and the Ki(1400) resonances, where the widths of the resonances are 
taken to be 90 and 174 MeV, respectively [jl4[. The shape of the r^— s>z/^ 7r^7r^7r+(7r°) 
background in the distributions is obtained from a data control sample. 
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Figure 7: (a) is the J^P distribution and (b), (c), and (d) are the invariant mass 
distributions of the data T~—>-i'rK~7i~K~^ candidate sample (points). The histogram 
represents the predicted distribution, in which the normahsation of the background and 
signal components comes from the results of the binned maximum likelihood fit to the 
four distributions. The shape of the r~— > z/,- K~7r~K"*" signal distribution is simulated 
in the fit by Monte Carlo events generated by the default version of Tauola 2.4. The 
shape of the r~— 7r~7r~7r"^(7r°) background in the distributions is obtained from a 
data control sample. 
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